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^ : ABSTRACT 

> : 

We discuss the influence of tfie ^^C{a,'y)^^0 reaction rate on the central He burning of stars in the mass 
p^ ' range 0.8 — 25 M0, as well as its effects on the explosive yields of a 25 Mq star of solar chemical composition. 

?H ' We find that the central He burning is only marginally affected by a change in this cross section within the 

currently accepted uncertainty range. The only (important) quantity which varies significantly is the amount 
of C left by the He burning. Since the ^■^C(o,7)^^0 is efficient in a convective core, we have also analyzed the 
influence of the convective mixing in determining the final C abundance left by the central He burning. Our 
main finding is that the adopted mixing scheme does not influence the final C abundance provided the outer 
border of the convective core remains essentially fixed (in mass) when the central He abundance drops below 
~ 0.1 dex by mass fraction; vice versa, even a slight shift (in mass) of the border of the convective core during 
the last part of the central He burning could appreciably alter the final C abundance. Hence, we stress that it 
is wiser to discuss the advanced evolutionary phases as a function of the C abundance left by the He burning 
rather than as a function of the efficiency of the ^^C(a, 7)^®0 reaction rate. Only a better knowledge of this 
cross section and/or the physics of the convective motions could help in removing the degeneracy between these 
two components. We also prolonged the evolution of the two 25 Mq stellar models up to the core collapse and 
computed the final explosive yields. Our main results are that the intermediate-light elements, Ne, Na, Mg and 
Al (which are produced in the C convective shell) scale directly with the C abundance left by the He burning 
because they depend directly on the amount of available fuel (i.e C and or Ne). All the elements whose final 
yields are produced by any of the four explosive burnings (complete explosive Si burning, incomplete explosive 
Si burning, explosive O burning and explosive Ne burning) scale inversely with the C abundance left by the He 
burning because the mass-radius relation in the deep interior of a star steepens as the C abundance reduces. 
We confirm previous findings according to which a low C abundance (~ 0.2 dex by mass fraction) is required 
to obtain yields with a scaled solar distribution. 
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1. Introduction 

"The rate of the ^^C(a,7)'^^0 during hydrostatic helium burning is of vital interest for 
explosive nucleosynthesis. It is this process that determines the abundances of ^^C and ^^O 
in the star, and thereby sets the stage for explosive burning... The rate is determined by the 
7.115 MeV level in the ^^O compound nucleus. At present the reduced width 6'^ of this 
resonance for a captures is not known." These sentences are taken from the 1973 issue of the 
Ann.Rev.Astron.Astrophys. and were written by Arnett to emphasize both the importance 
of this reaction in determining the final yields produced by the explosion of a supernova event 
and the fact that this rate was very uncertain. The experimental and theoretical efforts in 
the following 30 years led to constrain the reduced width 6^ of the 7.115 MeV level in ^^O 
and therefore the El component of the ^^C(q;,7)^^0 cross section. These studies, however, 
also pointed out how other components (not equally well constrained) contribute to the total 
cross section, so that the present determination of the stellar rate of the ^^C(q;, 7)^^0 is still 
affected by a large error. 

From an experimental point of view, in spite of the enormous efforts devoted to the 
measurement of this cross section, the corresponding rate at astrophysical energies is still 
far from being well established. The cross section around the Gamow peak is dominated by 
ground state transitions through four different processes: the two El amplitudes due to the 
low-energy tail of the 1~ resonance at Ecm = 2.42 MeV and to the subthreshold resonance 
at —45 keV, the E2 amplitude due to the 2+ subthreshold resonance at —245 keV and the 
direct capture to the ^^O ground state (plus the relevant interference terms). Besides ground 
state transitions, also cascades, mainly through the E2 direct capture to the 6.05 MeV 0"'" 
and 6.92 MeV 2^ states, have to be considered. Although they are believed to give a minor 
contribution (about 10 %) to the total cross section, no experimental data are available for 
such transitions. In the past twenty-five years many experiments have been set-up, most of 
them based on the detection of 7-rays produced by a captures in direct or inverse kinematics 
(Dyer & Barnes 1974; Redder et al 1985, 1987; Kremer et al 1988; Ouellet et al 1992, 1996; 
Roters et al 1999; Gialanella 2000). All these measures extend to a minimum energy of 
about 1 MeV and show systematic differences; below this energy, the extremely small value 
of the cross section (< 10 pb ) hampers the direct detection of 7-rays and extrapolation 
procedures have to be used in order to extract the astrophysical S-factor at the relevant 
energies (Eg = 300 keV for Tg = 0.18). Such an extrapolation, which is based on the fitting 
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of differential cross sections in tlie investigated region, requires also the inclusion of the phase 
correlation between the two incoming partial waves which contribute to the two multipoles. 

Additional information is provided by the elastic scattering data (Plaga et al 1987) and 
by the /5-delayed a-decay of ^^A^ (Buchmann et al 1993; Azuma et al 1994). Also the decay 
to the first excited state has to be included together to a possible non radiative EO ground 
state transition. As far as a consistent description of the El interference terms is concerned, 
it should be noted that the evaluation of the contribution of higher energy 1~ levels requires 
data to be taken at energies well above this resonance, where the competition with the 
background arising from the ^^C{a, n) reaction (or other neutron-producing reactions in 
inverse kinematics studies) makes cross section measurements very difficult. 

The above arguments make the extrapolated values of S(300) very uncertain. A global 
analysis (Buchmann et al 1996) of all the available data (surface fit) including the 7 decay 
which follows an a capture from ^^C, elastic scattering of a particles from ^^C and the a 
emission which follows a f3~ decay of ^^N (Buchmann et al 1993; Azuma et al 1994) yielded 
a wide range of results (from 62 keV b to 270 keV b ) for the extrapolated S-factor. The 
minimum and maximum values which bracket such a spread correspond to reaction rates 
(for Tg = 0.18) of 0.5-10"^^ and 2.2-10"^^ cm^/(mol-s), which can be compared to the data 
reported in the compilations of Caughlan and Fowler, 1988, (NaO"v = 0.8-10"^^ cm'^/(mol-s)), 
hereinafter CF88, and Caughlan et al., 1985, (Na(7v = 1.9 ■ 10~^^ cm^/(mol ■ s)), hereinafter 
CF85, which are generally used in stellar evolution calculations. Finally, a recent compilation 
(Angulo et al 1999) yields 0.9 ■ 10~^^ cm^/(mol • s) and 2.1 ■ 10~^^ cm^/(mol ■ s) as lower and 
upper bounds for this reaction rate, and 1.5 • 10^^^ cm^/(mol • s) as the recommended value. 

On the theoretical side, Arnett (1971) was the first to point out that the observed solar 
abundances of C and O could be used to limit the rate of the ^^C(q;,7)^^0. On the same 
line. Weaver and Woosley (1993) also tried to fix this rate by requiring the final explosive 
yields to have a scaled solar relative distribution. 

In addition to these efforts made to constrain this rate on the basis of the yields pro- 
duced, the direct infiuence of this process on the central He burning phase itself was also 
tested: in particular, Iben (1968, 1972) and Brunish & Becker (1990), by analyzing the behav- 
ior of a set of intermediate mass stellar models, found out that a change in the ^^C{a, 7)^^0 
reaction rate led to a change in the properties of the stars in the blue loop phase and hence, 
in turn, that it could modify the mass range capable of entering the Cepheids instability 
strip. Contrarily to these results, Umeda et al 1999, Zoccali et al (2000) and Bono et al 
(2000) found that a change in the ^^C(q;,7)^^0 rate does not modify the path of a star in 
the HR diagram. 
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For sake of completeness let us remind that also the properties of the cooling sequences 
of the White Dwarfs have been studied as a function of the relative abundances of C and O 
in the He exhausted core. We refer the reader to the papers by, e.g., Segretain et al (1994), 
Salaris et al (1997), Brocato et al (1999) and Chabrier et al (2000) for an overview of the 
main findings in this field. 

Though the partial effects of a change in this cross section on the evolution of a star 
have been addressed in several papers over the years (as we have already pointed out), 
a comprehensive and homogeneous analysis of its effects over an extended mass interval 
is still missing. Moreover, we believe that the interplay between the convection and the 
^^C{a, 7)^^0 in determining the chemical composition left by the He burning needs a deeper 
analysis. In this paper we will analyze the dependence of the central He burning phase on 
the ^^C(q;,7)^^0 reaction rate over an extended mass interval, together with its interplay 
with the convective mixing. We will also discuss the dependence of the final explosive yields 
produced by a 25 Mq on the C abundance left by the He burning. 

The paper is organized as follows: in the next section we briefly remind the main 
properties of the evolutionary code (FRANEC) adopted to perform all the computations. 
Section three is devoted to the discussion of the central He burnings of all our models while 
the advanced evolution of the 25 M0 stellar models are addressed in section four. The final 
explosive yields are presented in section five. A final discussion and conclusion will follow. 



2. The evolutionary Code 

All the evolutionary tracks have been computed with the latest release (4.8) of the 
FRANEC (Frascati RAphson Newton Evolutionary Code) whose earliest and latest version 
have been presented by Chieffi and Straniero (1989) and Chieffi, Limongi and Straniero 
(1998). All the latest available input physics have been adopted as discussed in Straniero, 
Chieffi and Limongi (1997). No mass loss has been taken into account. The network adopted 
in the present set of models includes 19 isotopes for the evolution of the low and intermediate 
mass stars and 179 isotopes for the evolutions of the 25 M© stars. Since we will discuss the 
effects of the overshooting and semiconvection on the stellar models, and since for historical 
reasons these words have been used to mean very different phenomena in stars of different 
mass, we brieffy remind what they refer to in the various mass ranges. 
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2.1. Overshooting and semiconvection in low mass stars 

During the central He burning, He is converted in C first and in O later. The increase 
of the C and O abundances in the convective core raises the opacity so that a jump in 
the radiative gradient forms at the border of the convective core. This is a condition of 
unstable equilibrium in the sense that the possible mixing (by whichever phenomenon) of 
the radiative layers just outside the border of the convective core would switch them from a 
stable to an unstable condition. The reason is that the C brought in the radiative layer raises 
the radiative gradient (through the opacity) so that it becomes intrinsically convective. This 
phenomenon, usually called "induced" overshooting, does not contain any free parameter 
which may be adjusted by hand since the process of "growth" of the convective core is fully 
controlled by the requirement that the positive difference between the radiative and adiabatic 
gradient cancels out. The word "induced" refers to the fact that this phenomenon is induced 
by the conversion of He in C and O. When the central He abundance drops below ~ 0.6 dex 
by mass fraction, the radiative gradient does not decrease any more monotonically moving 
outward but it forms a minimum well inside the formal border of the convective core. This 
occurrence triggers the formation of a region (outside the mass location corresponding to 
this minimum) in which the matter is only partially mixed: the condition which controls the 
degree of mixing occurring in this region is that the radiative gradient equals the adiabatic 
one (this equality is controlled, once again, by the opacity which, in turn, depends on the 
local abundances of C and O in these layers). This is the so called " semiconvective" region 
which forms in low mass stars. For a much more detailed discussion of these phenomena 
we refer the reader to, e.g., Castellani et al (1985). Since the "induced" overshooting and 
semiconvection completely depend on the fact that the opacity is strongly dependent on 
the chemical composition, it is clear that they become progressively less important, and 
eventually disappear, as the initial mass of a star increases because the electron scattering 
(which does not depend on the chemical composition in an environment deprived of H) 
becomes the main source of the opacity. In practice the semiconvective layer disappears for 
masses above ~ 5 Mq while the "induced" overshooting remains at least partially efficient 
up to ~ 20 Mq. 

The "real" existence of these phenomena in low mass stars is mainly supported by the 
star countings in the galactic globular clusters: in particular the ratio between the He burning 
stars (HB-stars) and those ascending the Giant Branch (the first and/or the second time) 
can be explained only if the central He burning timescale is the one obtained by including 
these two phenomena. Also in this case we refer the reader to Castellani et al (1985) for a 
careful discussion of these problems. 

During the latest part of the central He burning (i.e. when the He drops below 0.1 
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dex by mass fraction), it has been recognized that a runaway of the outer border of the 
convective core occurs (usually called Breathing Pulse -hereinafter BP- , see Castellani et al, 
1985 and Caputo et al, 1989): its main effect is that of engulfing fresh He towards the center 
and hence prolonging the central He burning lifetime. A discussion on the real existence of 
these instabilities is far beyond the purposes of this paper (but see Castellani at al. 1985); 
we simply want to stress the fact that their inclusion or suppression significantly alters also 
the abundances of C and O at the end of the He burning. 



2.2. Overshooting and semiconvection in massive stars 

The word overshooting is used, in this case, to mean the phenomenon which would 
allow the convective bubbles to penetrate the radiative layers surrounding a convective zone 
and hence to induce the mixing of a region larger (in mass) than classically allowed by the 
strict adoption of the Schwarzschild criterion. This is a mechanical phenomenon which is 
not confined to a specific evolutionary phase but which may be present at the border of any 
convective region. The extension of this overshooted region is, in principle, totally arbitrary 
and usually parameterized by imposing that the convective bubbles may reach a maximum 
extension over the formal convective border that is proportional to the pressure scale height 
(Hp). The existence of a convective core larger than permitted by the Schwarzschild criterion 
was invoked in the past in order to explain some observational data (see, e.g. Langer and 
Maeder, 1995, but see also Testa et al, 1999). Though we do not intend to discuss here the 
possible existence or not of a mechanical overshooting, it must be said that during the years 
the accepted size of this phenomenon in the central H burning phase progressively reduced 
from ^ 1 Hp down to less than 0.2 Hp. 

The word semiconvection is used, in this framework, to mean the partial mixing which 
(would) occur at the end of central H burning in the region of variable chemical composition 
left by the receding H-convective core in stars more massive than ~ 15Mq. When the star 
exhausts the H in the center and readjusts on a structure supported by an H burning shell, 
the radiative gradient overcomes the adiabatic one within these layers showing a gradient 
of chemical composition. While these layers would be definitely convective unstable if the 
Schwarzschild criterion were adopted to asses their stability, the adoption of the Ledoux 
criterion would maintain these layers stable. Observational constraints (see e.g. Langer 
and Maeder, 1995), mainly related to the observed number ratio between red and blue 
supergiants, seem to favor the Ledoux criterion, i.e. a partial or even negligible amount of 
mixing. 

Before closing this section let us clearly state that our standard computations are ob- 
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tained by adopting always the Schwarzschild criterion but in the central He burning phase 
where both the "induced" overshooting and semiconvection are properly taken into account 
while the BPs are quenched by forcing the central He abundance to be a monotonic not 
increasing function of time. No mechanical overshooting has been included. Moreover, no 
mixing is allowed in the semiconvective H-rich layers (which corresponds to a strict applica- 
tion of the Ledoux criterium). 



3. The central Helium Burning phase 

We followed the evolution of stellar model having 2.5 < M/Mq < 25, Y=0.285 and 
Z=0.02 from the Main Sequence up to the central He exhaustion. We also followed the 
central He burning phase of a typical Globular Cluster Horizontal Branch (HB) star, i.e. a 
star with an He core mass of 0.485 M©, a total mass of 0.6 M©, an initial He abundance 
Y=0.23 and a metallicity Z=0.001. All these evolutions have been computed twice: firstly 
by adopting the ^^C(q;,7)^^0 reaction rate provided by CF88 and secondly by adopting the 
one provided by CF85. 

Table 3 summarizes, for each mass, the main evolutionary properties in rows 1 to 12 
(each couple of columns refers to the values obtained with the CF88 and the CF85 rates, 
respectively). In order, left to right we report: the central He burning lifetime, the C and 
O abundances left by the He burning, the time spent by each model in the blue loop (i.e. 
at Log(Tcff > 3.80)), the fraction of the He burning lifetime spent in the blue loop, the He 
core mass at the He ignition, the maximum size of the convective core and the final He core 
mass at the He exhaustion. 

Figures 1 to 3 graphically show the effect of a change in the ^^C(q;,7)^^0 reaction rate 
from the CF88 to the CF85 one on the central He burning phase. 

In particular: Figure 1 shows the path followed by a selected sample of stellar masses 
in the HR diagram (the solid and dashed lines refer, respectively, to models computed with 
the CF88 and CF85 rate). Figure 2 shows, as filled dots, the difference in the He burning 
lifetimes (in percentage) obtained for the two rates as a function of the initial mass; Figure 
3 shows, instead, the percentage of the He burning lifetime spent in the blue loop (the filled 
and open dots refer, respectively, to models computed with the CF88 and CF85 rate). All 
these three figures show that an uncertainty of the ^^C(q;, 7)^^0 within the quoted range does 
not dramatically alter the "observable" properties of a star in the central He burning phase. 
In particular, the path followed by these stars in the HR diagram is practically unaffected 
by such a change whereas both the total He burning lifetime and the time spent in the blue 
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Fig. 1. — Path followed by selected models in the HR diagram; the dotted and solid lines 
refer, respectively, to the CF85 and CF88 cases. 
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Table 1. Main properties of the central Helium burning phase. 
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Fig. 2. — Dependence of the helium burning lifetime on the ^^C(q;,7)^^0 reaction rate as a 
function of the initial mass; each dot corresponds to (t^^^^ — t§^^^)/ty^^^. 
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Fig. 3. — Dependence of the "blue loop lifetime" on the ^^C(q;,7)^^0 reaction rate as a 
function of the initial mass; each dot shows the percentage of the helium burning lifetime 
spent at log Teg > 3.8; the open and filled dots refer, respectively, to the CF85 and CF88 
cases. 
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loop change by 10% at most. It is worth noting that the most massive star which experiences 
a blue loop in the central He burning phase changes from the 14 to 15 Mq as a consequence 
of the quoted change in the ^^C(q;,7)^^0 rate. 

Let us turn now to the chemical composition left by the He burning; Figure 4 shows the 
amount of C left by the He burning as a function of the initial mass. 

The filled symbols always refer to computations performed by adopting the CF88 value 
while the open symbols always refer to models computed by adopting the CF85 value; the 
dots refer to our "standard" models. The first thing worth noting is that the two sets of 
models show essentially the same dependence of the final C abundance on the initial mass 
and hence they are more or less systematically shifted one with respect to the other by 
0.20:0.25 dex. The general trend is that the C abundance left by the He burning increases 
as the initial mass reduces; a maximum is then reached for a mass of the order of 5 Mq 
and then a drop occurs for smaller values of the mass; the HB star behaves almost like the 
2.5 Mq. It has to be noted that the maximum variation of the C abundance is of the order 
of 0.15 dex over the full mass range under exam, and that this variation even reduces to 
0.1 dex for the masses larger than ~ 8 Mq. The existence of a smooth monotonic trend for 
masses larger than 5 Mq can be understood by reminding that a smaller mass favors the C 
production rather than its destruction because the 3a reaction rate scales with the square 
of the density while the ^^C(a;,7)^^0 scales linearly with the density. The inversion of the 
trend for masses smaller than ~ 5 Mq is probably due to the fact that stars with a very 
small He core mass spend enough time in the last part of the central He burning so that the 
conversion of C in O is strongly favored. The HB star behaves like the 2.5 Mq because they 
have a similar He core mass. 

Since the process we are dealing with occurs in a convective environment it is important 
to verify if, and to what extent, the final C abundance depends on the adopted convective 
scenario. Let us start with the standard one, i.e. the case in which the stability is controlled 
by the Schwarzschild criterion. Figure 5 shows, as a solid line, the typical behavior of 
the convective core as a function of the central He abundance. This figure shows that the 
convective core grows during the first part of the central He burning phase, reaches an 
asymptotic value and then remains constant (in mass) until the possible occurrence of the 
BP (if M < 15 Mq) or till the end of the central He burning (masses above the 15 Mq never 
develop BPs). Since in our standard scenario the BP are quenched out. Figure 5 represents 
the qualitative behavior of the convective core of all the stars in the mass interval here 
studied. As we have already mentioned above, the final C abundance which is obtained by 
adopting these assumptions is shown as filled dots in Figure 4. 

A second set of models spanning essentially the same mass interval has been recomputed 
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Fig. 4. — Carbon abundance left by the helium burning; the open and filled symbols refer, 
respectively, to the CF85 and CF88 cases. The dots represent our standard computations; 
the triangles represent the runs obtained by including 1 Hp of mechanical overshooting; the 
squares refer to the runs obtained by including the breathing pulses; the hexagon refers to 
the run in which 0.1 Hp of mechanical overshooting is imposed when the central helium 
abundance drops below 0.075 dex by mass fraction. 
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by imposing a large overshooting (1 Hp) during the central He burning phase. For sake of 
clarity let us remind the reader that the adoption of a large amount of overshooting auto- 
matically cancels out the possible formation of a semiconvective region because it completely 
mixes the region where the partial mixing should occur. In spite of a much larger (mass) 
size of the convective core (dotted line in Figure 5), the C abundance left by the He burning 
(shown as filled triangles in Figure 4) closely resembles the one obtained in the standard 
case. The reason is that the overshooting increases the size of the convective core but it 
does not alter the behavior of the border of the convective core, which remains essentially 
constant in mass during the latest part of the He burning; hence the run of both the central 
temperature and density as a function of the He abundance (see Figure 6) do not change 
significantly, as well as the rate at which He is converted in C, and C in O. The only effect 
of the overshooting is, in this respect, to increase the He burning lifetime as a consequence 
of the increased amount of available fuel. 

This picture changes drastically if one allows the border of the convective core to grow 
in mass when the central He drops below ~ 0.1 dex by mass fraction. This possibility 
may "naturally" occur if, e.g., one did not artificially dump out the occurrence of the BPs. 
The ingestion of fresh He in an environment very ^^C rich would favor, in this case, the 
^^C(a,7)^^0 rather than the 3a's, so that the final C abundance would be much lower than 
in the previous two scenarios. Moreover, since the number and the strength of the BPs scale 
inversely with the initial mass, it is clear that the lower the mass the larger will be their 
influence on the evolution of the star. 

The same effect which is obtained by means of the BPs may be obviously obtained in all 
cases in which even a small amount of fresh He is allowed to enter the convective core towards 
the end of the central He burning. In order to further stress how delicate the dependence of 
the final C abundance is on the behavior of the border of the convective core in the latest 
phases of the He burning, we show as a filled hexagon in Figure 4, the C abundance left by 
the He burning of a 15 M© in which just 0.1 Hp of mechanical overshoot is imposed when 
the central He burning drops below 0.075 dex by mass fraction: in this case the final C 
abundance even resembles the value obtained by adopting the CF85 rate. 

Before closing this section let us remark that, since massive stars do not have BPs and 
since the size of the convective core does not alter the final C abundance at the end of the 
He burning, one could be tempted to conclude that the C abundance left by the He burning 
depends only on the adopted value of the ^^C(q;,7)^^0 rate; however, since we do not feel 
confident to state that current uncertainties in the treatment of the convective core of the 
massive stars are merely confined to the size of the convective region itself, we prefer to 
conclude that in all the mass interval under exam the final C abundance left by the He 
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burning depends on both the mixing scheme and the adopted C(a,7) O reaction rate 



4. The advanced evolutionary phases of a 25 M0 

In the previous section we have shown the direct influence of the ^^C(q;,7)^^0 process 
on the central He burning of stars in a wide mass interval together with its interplay with 
the treatment of the convective core. The next logical step would be to follow the further 
evolution of all these stars in order to determine the final impact of this process on stars of 
different masses. Such a big project goes beyond the purposes of the present paper: in this 
section we will concentrate on the further evolution of the 25 Mq (taken as representative of 
the massive stars) up to the final collapse and explosion. The C abundance left by the He 
burning is ~ 0.4 for the CF88 rate and ~ 0.2 for the CF85 one. Since all the evolutionary 
properties discussed below depend directly on the C abundance left by the He burning but 
not (necessarily) directly on the adopted value for the ^^C(q;,7)^^0 rate (see the previous 
section), we think that the two runs obtained by adopting the CF88 and the CF85 rates 
must be discussed in terms of the C abundance left by the He burning (Cini). For this reason 
in this and the following sections we will change terminology: the run computed by adopting 
the CF88 rate will be referred to as the C0.4 case, to underline that the results directly 
depend on a C abundance equal to 0.424. Analogously, the run obtained by adopting the 
CF85 rate will be referred to as the C0.2 case. 

The main evolutionary properties of these two evolutions are summarized in table 4 and 
in Figures 7, 8 and 9. Table 4 reports, for each central burning, its lifetime, the size of the 
convective core, the abundance of the most abundant elements produced in the burning as 
well as the data relative to the convective shell episodes, if present. 

Figure 7 shows the path followed by the two stars in the Log(tc) -^ Log(pc) plane while 
Figures 8 and 9 show the behavior of the convective regions as a function of time. These 
figures summarize the temporal evolution of the two stellar models up to the time of the 
core collapse while a snapshot of the final structure at the time of the explosion is shown 
in Figures 10 and 11: the first one shows the internal run of the most abundant elements 
while the second one shows the final mass-radius relation together with the final electron 
mole density Yg. 

We will not discuss in detail the properties of the various burnings and we refer the 
reader to, e.g., Chieffi, Limongi and Straniero (1998) and Limongi, Straniero and Chieffi 
(2000) for a detailed analysis of the advanced burnings; here we want simply to underline 
how the C abundance left by the He burning, i.e. Cini, infiuences the advanced burnings. 
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Table 2. Main stages of the two 25Mq stars. 







CF85 






CF88 








H 


Burning 




TH(.yr) 




5.81(6) 






5.81(6) 


Mcc(Mq) 




12.7 






12.7 








H. 


e Burning 




At(H-exh.He- 


■ign.) 


2.70(4) 






2.70(4) 


THe{yr) 




5.8(5) 






6.37(5) 


Mcc {Mq ) 




5.6 






5.8 


^'^He conv sh' 


eiiiyr) 


1.6(4) 






1.5(4) 


^Mue conv a 


helliMQ) 


2.1 






2.2 


12c 




0.424 






0.200 


Wq 




0.546 






0.769 


C Burning 


At(He-exh.C- 


■ign.) 


1.17(4) 






1.03(4) 


Tc(,yr) 




5.76(3) 






4.56(3) 


Mcc{Mq) 




0.5 








^tlC conv sh, 


5u(yr) 


91 






1 


^Mic conv s 


Mi (MO) 


1 






1.2 


^hc conv sh, 


.H(yr) 


40 






0.2 


AM2C conv s 


hcn{M&) 


3 






2.4 


16q 




0.378 






0.674 


20Are 




0.478 






0.260 


24Mg 




0.014 






0.076 


Ne Burning 


-rjvefor) 




37.9 






6.01 


Mcc{Mq) 




0.56 






0.77 


16(3 




0.632 






0.810 


24Ms 




0.139 






0.072 


28 Si 




0.143 






0.071 











Burning 




To(yr) 




1.62 






0.274 


Mcc{Mq) 




1.26 






0.98 


28 5i 




0.561 






0.604 


325 




0.014 






0.008 


34^ 




0.336 






0.150 








Si 


Burning 




TSiiyr) 




0.21 






0.0167 


Mcc{Mq) 




0.95 






1.28 


56 Fe 




0.507 






0.674 


eoAfj 




0.007 






0.023 
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Let us firstly note that the region outside the CO core, i.e. the He and H rich layers, is 
not significantly infiuenced by Cini because the typical timescale on which this outer region 
evolves is in any case much longer than the lifetime of all the advanced burning phases put 
together. 

The evolution of the CO core, on the contrary, will largely depend on Cjni since both its 
physical and chemical evolution will depend on the amount of fuel available in the C burning 
(both central and shell burnings). Since, as it is well known, the neutrino losses become a 
very efficient energy sink when the central temperature raises above ~ 8 x 10^ K, and since 
the formation of a convective core requires the nuclear energy (which depends quadratically 
on the C abundance) to overcome the neutrino losses, it is clear that a convective core may 
form in the central C burning phase only if Cini is larger than a threshold value. In our case 
a convective core forms in the C0.4 run while C burns in a radiative core in the C0.2 run. 
Once the C is exhausted in the center, the following C shell burning occurs (in both cases) 
through the formation of successive convective episodes. In spite of the very different amount 
of available fuel and of the details of the shell evolution, the last C convective shells obtained 
in the two cases show some conspicuous similarities: in particular the outer border of the 
convective shell is essentially insensible to Cini because it is fixed by the location of the He 
shell (which is located at the same mass coordinate in both cases) while the inner one is only 
mildly dependent on Cini in the sense that the location of the burning shell (which marks 
the base of the convective shell) is slightly shifted outward in the run with the lowest initial 
C abundance. Roughly speaking, the size of the convective shell reduces by almost 20% by 
mass fraction by increasing the initial C abundance from 0.2 up to 0.4 dex. The other very 
important similarity between the two runs is that, in spite of the very different amount of 
C present in the two convective shells, both models burn almost completely the C present 
in the shell. The existence of these similarities imply that the final chemical composition 
within the convective shell largely depends on Cini. The reason is obviously that, since the C 
is almost completely destroyed in both cases and since the mass size of the convective shell 
is similar, the abundances of the elements mainly produced by the C burning will directly 
depend on the available fuel, i.e. on Cini. 

The region behind the C burning shell continues to contract (and to heat) in order to 
counterbalance the energy losses and hence to further manipulate the chemical composition 
(through the Ne, O and Si burnings) up to the time of the collapse. In order to understand 
how the yields coming from this internal region depend on Cini it is not necessary to discuss 
in detail the various burnings beyond the C one but simply to understand how the final 
mass-radius relation depends on Cini. 

Figure 12 shows the mass-radius relation relative to the C0.2 (thin lines) and to the C0.4 



-17- 



(thick lines) cases at three selected points: the solid lines mark the end of the central C 
burning, the dotted lines refer to the beginning of the central Ne burning while the dashed 
ones refer to the last computed model. A comparison between the two structures shows that 
the two models reach the end of the central C burning with a similar M-R relation (within 
the first four solar masses). During the further evolutionary phases which eventually lead to 
the core collapse, the very interiors of the two stars (M < 1.4 Mq) continue to contract by 
maintaining a similar M-R relation, while the more external regions reach the time of the 
core collapse with very different M-R relations. The lines showing the two M-R relations 
at the beginning of the central Ne burning phase reveal that most of the difference actually 
start before the central Ne burning phase. In the time interval which elapses between the 
end of the central C burning and the beginning of the central Ne burning the CO core 
experiences a phase of contraction in which it gains from the gravitational field the amount 
of energy necessary to maintain the hydrostatic equilibrium. In this transient phase, the 
energy requirement by the CO cores is partially alleviated by the formation of one (or more) 
convective C shell episodes. These convective episodes stop for a while the advancing C 
burning front and allow the C-burning shell to burn a "reservoir" of fuel while remaining 
essentially fixed in mass: such an occurrence helps in slowing down the contraction of the 
region above the C burning front. The larger amount of C available in the C0.4 case allows 
a more effective support of the layers above the C burning front and hence the formation 
of a M-R relation less steeper than in the other case: the dotted lines in Figure 12 clearly 
show such an occurrence. Though both models will further strongly contract up to the time 
of the core collapse, the differences in the M-R relations which form before the Ne ignition 
remain till the final explosion. 



5. The explosive yields 

Once a pre supernova model is obtained, it is necessary to simulate in some way the 
explosion in order to compute the final yields. We refer the reader to the papers by Limongi, 
Straniero and Chiefii (2000) and by Limongi, Chiefii and Straniero (in preparation) for a 
comprehensive discussion of the technique we adopt to simulate the passage of a shock 
wave: here it suffices to say that the two explosions have been followed by assuming that 
a shock wave successfully escapes the iron core giving a final kinetics energy of 1.2 x 10^^ 
ergs and that the mass cut has been arbitrarily chosen to eject 0.05 Mq of ^^Ni. In order 
to understand the dependence of the explosive yields on the C abundance left by the He 
burning it is important to remind a few key properties of the explosion. Once the shock 
wave generated by the rebounce of the core escapes the iron core, it moves through the 
mantle of the star without loosing essentially any energy: hence the peak temperature at 
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the shock front lowers as the shock moves outward simply because it expands adiabatically 
and not because it crosses progressively larger layers of the star. This means that the peak 
temperature at the shock front is a function of its geometrical distance from the center and 
not of the amount of mass crossed by the shock wave. Hence, once the energy of the shock 
front exiting the iron core is fixed, it is possible to determine " a priori" the radii at which the 
various peak temperatures will be reached. By reminding that at a good approximation it 
can be assumed that E = 4/37rR^aT^, it can be easily determined that, for an initial energy 
of 1.2 X 10^^ ergs, a peak temperature of 5 x 10^ K is reached at r=3900 Km, a T = 4 x 10^ K is 
obtained at r=5300 Km, a T = 3.3 x 10^ K is obtained at r=6800 Km while T = 1.9 x 10^ K is 
reached at r= 14200 Km. This grid of radii corresponding to these key temperatures defines 
the volumes of space within which the matter will be exposed to, respectively, complete 
explosive Si burning (Tg > 5), incomplete explosive Si burning (5 > Tg > 4), explosive O 
burning (4 > Tg > 3.3) and C and Ne explosive burnings (3.3 > Tg > 1.9). Apart from the 
C and Ne explosive burnings, all the other three explosive burnings leave a specific (i.e. per 
unit mass) chemical composition which depends on the pre explosive chemical composition 
only through its local degree of neutronization (which may be expressed, e.g., by means of 
the electron mole density Yg). This means that a change in the C abundance left by the 
He burning does not modify the specific yields produced by these explosive burnings (the 
degree of neutronization reached by the matter in these zones is mainly determined by the 
conversion of ^^N - which means the initial Zcno - in ^^Ne). Hence Cini influences the final 
yields of the elements produced by these burning only through its influence on the final M-R 
relation (which means, in practice, the amount of matter located in the various key zones). 
Table 5 shows, for both runs, the amount of matter exposed to the three explosive burnings: 
In accordance with the M-R relations obtained in the two cases, the amount of matter 
exposed to both the explosive Oxygen burning and the incomplete explosive Si burning is 
significantly larger in the C0.2 case. Only the amount of matter exposed to the complete 
explosive Si burning is larger in the C0.4 case. This is, however, simply the consequence of 
the chosen mass cut: in fact, the final pre-explosive structure obtained in the C0.2 run is so 
compact that the required amount of ^^Ni is already almost completely synthesized by the 
incomplete explosive Si burning. 

Keeping in mind these properties of the explosion we can now turn to the analysis of 
the dependence of the explosive yields on Cini. 

Table 5 shows the isotopic yields (2.5 x 10'^ s after the rebounce) produced in the two 
cases while table 5 and Figure 13 show a comparison between the elemental (fully decayed) 
yields. In the following we will focus our attention only on the elemental yields. This will be 
equivalent, in general, to speak about the most abundant isotope: of course if more than one 
isotope is important to describe the behavior of an element we will address all the important 
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Fig. 5. — Template behavior of the border of the convective core versus the central He 
abundance for the present set of models. The solid line refers to a standard model while the 
dotted line refers to the same model but computed by adding 1 Hp of overshooting. 



Table 3. Mass intervals exposed to the various explosive burnings 



zone 


AM(Co.2) 


AM(Co.4) 




(Mo) 


(Mo) 


St. 





0.03 


Slix 


0.15 


0.09 


o. 


0.22 


0.13 


{CkNe)x 


0.82 


0.68 
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Fig. 6. — Run of the central temperature (upper panel) and density (lower panel) as a 
function of the central He abundance. The solid line refers to the standard 10 M©, while the 
dotted one refers to the 10 M© computed with 1 Hp of mechanical overshooting. The two 
dashed lines show, as a reference, the behavior of a 12 M© and an 8 Mq. All these models 
were computed by adopting the ^^C(q;,7)^^0 reaction rate as provided by CF88. 
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Fig. 7. — Run of the two 25 Mq in the Log(Tc) -r- Log(pc) plane. The sohd and dotted hnes 
refer, respectively, to the C0.4 and C0.2 cases. See text. 
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Fig. 8. — Temporal (properly adapted) behavior of the convective zones which form during 
the evolution of the 25 Mq stellar model in the C0.2 case. 
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Fig. 9. — Temporal (properly adapted) behavior of the convective zones which form during 
the evolution of the 25 Mq stellar model in the C0.4 case. 
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Fig. 10. — Structural profiles of the most abundant isotopes within the two test 25 M© 
stellar models at the time of the core collapse. 
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Fig. 11. — Comparison of the two final Mass- Radius and Ye- Radius relations. 
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Fig. 12. — Mass- Radius relation at selected points along the evolution: end of the central C 
burning (solid lines), beginning of the central Ne burning and (dotted lines) and last model 
(dashed lines). The thin lines refer to the C0.2 case while the thick ones refer to the C0.4 run. 
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Table 4. Isotopic yields 2.5 x 10^ s after the rebounce. 





Co 


.2 


Co.4 




Co 


2 


Co.4 




(M©) 


(Mo) 




(M©) 


(M©) 


H 


1.04 


X 


10+1 


1.04 


X 10+1 


39k 


2.07 


X 


10-4 


1.58 


X 10-4 


2H 


2.20 


X 


10-16 


2.19 


X 10-16 


40k 


3.04 


X 


10-6 


4.14 


X 10-6 


3He 


2.79 


X 


10-1 


2.79 


X lO-'' 


41k 


2.18 


X 


10-5 


1.78 


X 10-5 


"He 


8.02 


X 


lOO 


8.13 


X 10" 


40 Ca 


1.30 


X 


10-2 


6.83 


X 10-2 


6Li 


1.33 


X 


10-9 


1.33 


XlO-9 


42 Ca 


8.36 


X 


10-5 


5.16 


X 10-5 


^Li 


6.74 


X 


10-11 


6.79 


X 10-11 


43 Ca 


5.12 


X 


10-6 


5.94 


X 10-6 


9 Be 


3.39 


X 


10-10 


3.40 


X 10-1" 


44 Ca 


4.60 


X 


10-5 


5.76 


X 10-5 


10b 


2.25 


X 


10-9 


2.25 


X 10-9 


46 Ca 


1.74 


X 


10-6 


2.77 


X 10-7 


"B 


2.07 


X 


10-8 


2.08 


XlO-8 


48 Ca 


3.04 


X 


10-6 


3.11 


X 10-6 


12c 


5.04 


X 


10-1 


6.83 


X 10-1 


45 Sc 


4.51 


X 


10-6 


3.31 


XlO-6 


13c 


2.22 


X 


10-3 


2.24 


X 10-3 


46 'pj 


3.12 


X 


10-5 


2.13 


X 10-5 


14N 


7.87 


X 


10-2 


8.11 


XlO-2 


47rpj 


6.48 


X 


10-6 


6.83 


X 10-6 


15n 


2.70 


X 


10-5 


2.71 


X 10-5 


48 Ti 


1.65 


X 


10-4 


1.18 


X 10-4 


16o 


2.39 


X 


lOO 


1.71 


X 10" 


49^1^; 


1.74 


X 


10-5 


1.43 


X 10-5 


i^O 


1.29 


X 


10-1 


1.30 


X lO-'' 


5"Ti 


1.28 


X 


10-5 


1.45 


X 10-5 


18o 


5.04 


X 


10-4 


2.56 


X 10-"' 


50v 


1.48 


X 


10-7 


1.52 


X 10-7 


19p 


1.27 


X 


10-5 


1.09 


X 10-5 


Sly 


3.04 


X 


10-5 


2.23 


X 10-5 


20Ne 


3.53 


X 


10-1 


1.02 


X 10" 


50 Cr 


1.79 


X 


10-4 


1.25 


X 10-4 


21 Ne 


2.18 


X 


10-3 


1.83 


X 10-3 


52 Cr 


2.48 


X 


10-3 


1.44 


X 10-3 


22 Ne 


5.93 


X 


10-2 


5.30 


X 10-2 


53 Cr 


2.71 


X 


10-4 


1.78 


X 10-4 


23Na 


1.60 


X 


10-2 


3.23 


XlO-2 


54Cr 


3.24 


X 


10-5 


3.53 


X 10-5 


24Mg 


8.21 


X 


10-2 


2.99 


X 10-1 


55 Mn 


1.22 


X 


10-3 


8.79 


X 10-3 


25 Mg 


2.40 


X 


10-2 


3.41 


X 10-2 


54Fe 


1.41 


X 


10-2 


1.01 


X 10-2 


26 Mg 


1.86 


X 


10-2 


2.61 


XlO-2 


56Fe 


7.36 


X 


10-2 


7.38 


XlO-2 


27A1 


1.26 


X 


10-2 


3.20 


XlO-2 


57Fe 


1.74 


X 


10-3 


2.28 


X 10-3 


28 Si 


2.09 


X 


10-1 


1.87 


X 10-1 


58Fe 


1.06 


X 


10-3 


1.15 


X 10-3 


29 g; 


6.06 


X 


10-3 


7.65 


X 10-3 


59 Co 


4.49 


X 


10-4 


6.48 


X 10-4 


30 g; 


6.69 


X 


10-3 


7.57 


X 10-3 


58Ni 


2.17 


X 


10-3 


3.32 


X 10-3 


31p 


1.78 


X 


10-3 


2.20 


X 10-3 


60Ni 


9.36 


X 


10-4 


1.02 


X 10-3 


32 g 


1.06 


X 


10-1 


7.00 


X 10-2 


6iNi 


2.04 


X 


10-4 


2.30 


X 10-4 


33 g 


5.72 


X 


10-1 


6.05 


X 10-'' 


62Ni 


5.50 


X 


10-4 


6.62 


X 10-4 


34g 


6.74 


X 


10-3 


7.06 


X 10-3 


64Ni 


4.96 


X 


10-4 


4.83 


X 10-4 


36 g 


2.83 


X 


10-5 


2.99 


X 10-5 


63 Cu 


2.85 


X 


10-4 


3.04 


X 10-4 


35 ci 


2.45 


X 


10-1 


2.43 


X 10-'' 


65 Cu 


1.20 


X 


10-4 


1.40 


X 10-4 


37C1 


1.92 


X 


lo-* 


2.10 


X 10-" 


64Zn 


5.61 


X 


10-5 


1.14 


X 10-4 


36 Ar 


1.64 


X 


10-2 


9.40 


X 10-3 


66 Zn 


1.71 


X 


10-4 


2.65 


X 10-4 


38 Ar 


3.03 


X 


10-3 


2.25 


X 10-3 


67Zn 


2.98 


X 


10-5 


5.13 


X 10-5 


40 Ar 


6.37 


X 


10-6 


6.06 


X 10-6 


68 Zn 


4.97 


X 


10-4 


7.04 


X 10-4 
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Table 5. Elemental yields fully decayed. 





Co. 2 


Co. 4 




(Mq) 


(Mq) 


H 


1.04 


X 10+01 


1.04 


X 10+01 


He 


8.14 


X 10+00 


8.02 


X 10+00 


C 


6.85 


X lO-oi 


5.06 


X 10-01 


N 


8.11 


X 10-02 


7.88 


X 10-02 


O 


1.71 


X 10+00 
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Fig. 13. — The logarithmic ratio between the Yields produced by the C0.4 run and those 
produced by the C0.2 run. 
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ones. A proper discussion of Figure 13 requires the knowledge of the production site of 
each element. Schematically, four main groups of elements may be identified (the isotopes 
within the brackets at the right side of each element show the main isotopes, at least in these 
runs, which determine the final elemental yields): the first one, which includes Ne (^°Ne), 
Na (23Na), Mg (^^Mg), Al (^^Al), P (^ip), CI (^^Cl and ^^Cl) and Sc (^^Sc and ^^Ca), is 
produced in the C convective shell; the second one (i.e. the "golden" group, see Limongi, 
Chiefii and Straniero in preparation) is produced by both the incomplete explosive Si burning 
and the explosive O burning and includes Si (^^Si), S (^^S), Ar (^^Ar), Ca (^"^Ca) plus K 
(^^K) which is synthesized only by the explosive O burning); the third one is produced only 
by the explosive incomplete Si burning and includes Ti (^^Cr), V (^^Cr), Cr (^^Cr, ^^Mn and 
^^Fe) and Mn (^^Mn, ^^Fe and ^^Co); the last one is produced by the complete explosive Si 
burning and includes Fe (^^Ni, ^^Fe and ^^Fe), Co (^'^Co), and Ni (^^Ni, ^^Ni, ^^Ni, ^^Ni and 
^^Ni) (We can't discuss Cu and Zn because they are just the upper end of our network). Iron 
is actually produced also as ^^Ni by the incomplete explosive Si burning. The light elements 
H to F will be discussed separately. Note that all the following discussion strictly holds only 
for the 25 Mq even if it probably may be considered more or less valid in general. 

Four out of the seven elements pertaining to the first group, namely Ne, Na, Mg and 
Al, present a very similar behavior: they are produced in the C convective shell and are 
partially destroyed by the explosion. Hence their final yields largely depend on the pre 
explosive ones. Figure 14 shows, as a typical example, the Ne profile at the time of the core 
collapse (dashed line) and after the passage of the shock wave (solid line). All four elements 
show a similar dependence on Cini, in the sense that they scale more or less uniformly 
(and directly) with the C abundance left by the He burning (see Figure 13). However this 
occurrence is somewhat accidental. Ne and Mg are direct products of C burning and hence 
it is quite reasonable that they scale similarly with Cjni. Na, on the contrary, though it is 
a primary product of C burning, settles rapidly at its equilibrium value between production 
and destruction: this equilibrium value is almost independent on the initial C abundance 
while it largely depends on the temperature in the sense that its abundance increases as the 
temperature decreases. Since the temperature at the base of the C convective shell scales 
inversely with Cini if follows that the Na yield increases as Cini increases (hence behaving 
similarly to Ne and Mg). What it is quite accidental is that it scales even quantitatively in 
a way similar to that of Ne and Mg. Al shows an even different behavior: it has a secondary 
origin (i.e. it descends mainly from the initial abundance of CNO) and it is formed mainly 
by the sequence ^^Ne(Q;, n)^^Mg(n, 7)^^Mg(p,7)^'^Al plus the additional (primary) sequence 
^^Na(Q;, p)^^Mg(p, 7)^'^A1. Since the first (dominant) sequence originates from the ^^Ne (plus 
^^Mg and ^^Mg), one could expect the Al yield to be controlled first of all by the amount of 
^^Ne +^^ Mg +^^ Mg, i.e. the seed nuclei. This is not the case because the starting reaction 



-30- 



of this sequence requires the presence of a particles (which in this conditions come directly 
from the ^^C +^^ C reactions); since, even for the largest initial C abundances, not all the 
^^Ne +^^ Mg +^^ Mg is fully converted in Al, it happens that the final Al production scales 
with the a's available, which means that it scales with the initial C abundance. Note that, 
as the initial metallicity reduces, the second (primary) channel becomes more important 
so that the final Al abundance always depends on Cinj and hence it behaves as a primary 
element. Let us now turn to the other three elements pertaining to the first group: i.e. P, 
CI and Sc. Each of them has a specific history and hence they will be discussed individually. 

P is produced by the Ne burning and its site of hydrostatic production is deep enough 
that it is completely destroyed by the passage of the shock wave: the P which will be ejected 
by the explosion is only the one directly produced by the shock itself when passing through 
the Ne rich layers: its typical final profile is shown in Figure 15. It is possible to identify 
two main chains of processes which lead to the production of P, the first one, primary, starts 
directly from the ^'^Ne and leads through several sub channels to P, and a second one which 
starts from ^^Mg and hence it is a typical secondary sequence. The not negligible presence 
of a true secondary component which does not depend on Cini leads to a rather mild global 
dependence of P on Cini. 

CI is a complex element since both stable isotopes, ^^Cl and ^^Cl, contribute significantly 
to the CI elemental yield. ^^Cl is similar to P because it is produced by the explosive burning 
while ^^Cl is produced in the central He burning phase and it is then preserved down to the 
base of the C convective shell up to the time of the core collapse. The passage of the shock 
wave partially destroys the ^^Cl, in a manner similar to Ne. Since none of the two isotopes 
descends directly from Cini, the final abundance of this element is practically unaffected by 
a change in Cini. 

Sc, the last element pertaining to this group, scales inversely with Cini. Its production 
depends on the abundance of ^^Sc itself and *^^Ca which decays in ^^Sc. Such a scaling 
may be explained by noting that *^^Ca is a branching point in which two concurrent process 
compete, i.e. the /9 decay and the neutron capture, and hence that the final fate of the ^^Ca 
depends on the neutron density; this last quantity scales directly with the C abundance 
simply because the same total amount of neutrons is released on different timescales. In the 
Co. 2 case the neutron flux is large enough that most of the matter coming from *^^Ca goes 
to *^^Ca first and to "^^Ca later: the abundance of *^^Ca is in this case high and determined 
by the equilibrium condition between ^^Ca(n, 7) and '^^Ca(n, 7). In the C0.4 case, on the 
contrary, the low neutron flux allows the ^^Ca to decay in '^^Sc, so that the final abundance 
of ^^Sc is settled by the competition between '^^Ca(n, 7) and '^^Sc(n, 7). 

Four out of the five elements pertaining to the Golden group, namely Si, S, Ar and 
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M(MJ 

Fig. 14. — Neon profile (C0.4 case) within the star at the time of the explosion (dashed line) 
and after the passage of the shock wave (solid line). 
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Fig. 15. — Phosphorus profile (C0.4 case) within the star at the time of the explosion (dashed 
line) and after the passage of the shock wave (solid line). 
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Ca, have a similar history. Since they are synthesized by both the incomplete explosive 
Si burning and the explosive O burning, their yields do not depend on the pre explosive 
composition (apart from the neutron excess which is however the same in the two runs) but 
only on the amount of mass exposed to these burnings. Since, as we have shown above, the 
lower the C abundance the larger is the amount of mass exposed to these explosive burnings, 
the final yields of these elements scale inversely with Cini. The relative abundances of these 
elements change mildly, but systematically, by a change in Cini so that the four elements (Si, 
S, Ar and Ca) are progressively more overproduced going from the C0.4 to the C0.2 case. 

The yields of the elements produced by the incomplete explosive Si burning depend on 
the amount of mass exposed to a peak temperature in the range 4-3 billions degree and hence 
they scale inversely with Cini. The opposite behavior of Ni and Co can be easily understood 
by noting that they are mainly produced by the complete explosive Si burning and that 
essentially no matter coming from this region is ejected in the C0.2 case (this is obviously 
true only for this specific choice of the mass cut). 

Let us now address also the dependence of the yields of the three long-lived radioactive 
isotopes, ^^Al, ^*^Fe and ^^Ti on Cini because of their importance as gamma ray emitters. 
These three isotopes are produced by the explosion itself, since their pre explosive abundance 
is either negligible or destroyed by the passage of the shock wave. The first two, i.e. ^^Al 
and ®°Fe, are produced in the region where the peak temperature of the shocks reaches a 
value of the order of 2 x 10^ K. The amount of ^^Al produced depends on the abundances 
of ^^Mg, ^^Na and ^°Ne (or ^^C) because the a particles produced by the Ne (and/or C) 
burning are captured by ^^Na which releases the protons which are then partially captured 
by ^^Mg to produce ^^Al. Hence the yield of this isotope scales directly with Cini and varies 
from 3.28 x 10"^ (for C0.4) to 2.67 x 10"^ (for Co.2)- The production of ®°Fe requires a 
double neutron capture on the stable isotope ^^Fe. Neutrons are mainly produced by the 
capture of a particles on ^^Ne and hence it is the abundance of this isotope at the time of 
the core collapse (in the region where the peak temperature is of the order of 2 x 10^ K) 
which controls the yield of the ^"^Fe. Since the final abundance of ^^Ne scales inversely with 
Cini, also the final ^'^Fe abundance will follow the same trend. In particular we predict an 
^°Fe yield equal to 4.94 x lO'^ (C0.4 case) and to 2.92 x IQ-^ (C0.2 case). ^^Ti is produced by 
the complete explosive Si burning and hence its yield will depend, among the other things, 
on the mass cut and on the degree of freeze-out experienced by the most internal layers of 
the ejecta. Under the (arbitrary) assumption that in both runs 0.05 M0 of ^^Ni are ejected, 
we obtain that the yield of this isotope scales inversely with Cini, and in particular that it 
reduces from 1.1 x 10"^ to 3.4 x 10"^. 

The light elements C and O are only marginally affected by the explosion and hence 
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they reflect essentially their pre explosive abundances: it goes without saying that while the 
C yield scales directly with Cini, the O yield scales inversely with Cinj. N is produced in the 
H burning and hence it is not affected by Cini • Fluorine is synthesized in the He burning shell 
by the chain -^^N(q;,7)^^F(/3)-'^^0(p, q;)^^N(q;,7)^^F and scales inversely with Qni because of 
both the mass size of the He shell and the ^^O abundance scale inversely with Cini- 



6. Final discussion and conclusions 

In the previous section we have shown that the elements between Ne and Ni may be 
divided in groups of nuclei which behave more or less similarly with respect to a change in 
Cini- Roughly speaking we can say that while the elements produced in the C convective 
shell scale directly with Cini (the larger the C abundance left by the He burning the larger 
the final abundance of these elements), most of the elements synthesized by the explosive 
burnings scale inversely with Cini because of the steeper M-R relation. Hence a comparison 
between the solar chemical composition and the yields obtained with the two values of Cini 
could help to constrain the real abundance of C left by the He burning. However, in order 
to obtain a robust comparison it would be necessary to integrate at least over a stellar 
generation extending between 13 and 30 Mq; at present we do not have such an extended 
set of computations for two different values of Cini. Nonetheless we think that it is in 
any case interesting to show a simple comparison between our data and the solar chemical 
distribution because, as it has already been noted several times (e.g. Woosley and Weaver 
(1982)), the star which influences more pronouncedly the chemical composition of the ejecta 
of a generation of stars is of the order of the 25 Mq if the adopted IMF is the Salpeter one. 
Figure 16 shows the production factors of all the elements discussed above for the two values 
of Cini: the filled squares refer to the C0.4 case while the open dots refer to the C0.2 one. 

By the way, the production factor is defined as the ratio between the amount of mass 
(in solar masses) ejected as a given element and the amount of mass (in solar masses) the 
same elements would have if all the ejecta would have a solar chemical composition. It goes 
without saying that all the elements sharing a similar production factor maintain scaled solar 
relative proportions. 

In the C0.4 case it can be seen that, even if slightly, C is overproduced with respect to 
Oxygen; this would imply that, at the very least, there would not be room for C production 
by other kinds of stars. Elements Ne to Ca show a production factor which systematically 
reduces with the atomic number Z with also large deviations from the O level: for example 
the block of elements Ne to Al are overproduced by a factor of 3-5 with respect to O. If this 
were correct, it would automatically mean that O would not be mainly produced by massive 
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stars (obviously if one assumes the elements from Ne to Al to be produced by massive stars). 
The elements from Sc to Ni are all under-overabundant with respect to O (but the Co): this 
result is in line with the idea that these elements probably come from the ejecta of a Type 
la Supernova. The only exception is Co which has a [Co/0] > 0: if this were correct we 
would face the unpalatable situation that the bulk of the Iron peak nuclei would come from 
one kind of star (type la Supernova) while just one single element of this group, Co, would 
come from massive stars. 

In the Co.2 case the only element which is clearly a problem is Na, which is embarrass- 
ingly overproduced with respect to O ([Na/0] ~ 0.3). All the other elements are more or 
less lined up with the currently mostly accepted scenario: all the elements from O to Ca 
share a very similar production factor, which means that they all come from massive stars. 
C and N are underproduced with respect to O so that other sources (AGB's and the like) 
may contribute to their synthesis. The same occurs for the elements beyond Ca which are 
all underproduced with respect to Oxygen by a factor 3-4, leaving wide room for a Type la 
contribution to the galactic enrichment. 

We hence conclude that, if the 25 Mq may be considered the leading polluter of the 
interstellar medium and if the solar chemical composition is the "reference" distribution, 
a low C abundance, of the order of 0.2 dex by mass fraction, should be left by the He 
burning. A result very similar to the present one was already obtained by Weaver and 
Woosley (1993), who computed a large set of evolutionary models over different values of the 
^^C(q;,7)^^0: by comparing their results to the solar distribution they concluded that the 
"correct" ^^C(q;, 7)^^0 rate should be of the order of 1.7 times the rate quoted by CF88. Note 
that the C abundance they obtain at the end of the He burning in the 25 M© by adopting 
their "best" rate is 0.18 dex, i.e. remarkably similar to our standard case. However, we have 
shown that the final C abundance left by the He burning does not depend on the size of the 
convective core only if its border remains constant in mass when the central He drops below 
~ 0.1 dex; if this were not the case, the final C abundance would strongly depend on the 
behavior of the convective core. Since we do not feel confident to state that we can robustly 
model the behavior of the convective core, we prefer to interpret both our and their results 
in terms of the C abundance left by the He burning rather than in terms of an effective 
^2C(a, 7)^60 rate. 

More light will be shed shortly on this topic by an European task force which is beginning 
a new measurement of such a tough cross section. Stay tuned... 
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Fig. 16. — Comparison between the production factors obtained in the C0.4 case (filled 
squares) and those produced in the C0.2 case (open dots). 



